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Context. Star formation theories are currently divergent regarding the fundamental physical processes that dominate the substellar 
regime. Observations of nearby young open clusters allow the brown dwarf (BD) population to be characterised down to the planetary 
mass regime, which ultimately must be accommodated by a successful theory. 

Aims. We hope to uncover the low-mass population of the p Ophiuchi molecular cloud and investigate the properties of the newly 
found brown dwarfs. 

Methods. We use near-IR deep images (reaching completeness limits of approximately 20.5 mag in 7, and 18.9 mag in H and K^) 
taken with the Wide Field IR Camera (WlRCam) at the Canada France Hawaii Telescope (CFHT) to identify candidate members of 
p Oph in the substellar regime. A spectroscopic follow-up of a small sample of the candidates allows us to assess their spectral type, 
and subsequently their temperature and membership. 

Results. We select 110 candidate members of the p Ophiuchi molecular cloud, from which 80 have not previously been associated 
with the cloud. We observed a small sample of these and spectroscopically confirm six new brown dwarfs with spectral types ranging 
from M6.5 to M8.25. 

Key words, stars: formation ~ stars: low-mass, brown-dwarf - stars: planetary system 



1. Introduction 

The determination of the initial mass function (IMF) across 
the entire stellar and substellar mass spectrum is a funda- 
mental constraint for star formation theories (see, for exam- 
ple, Bonnell et al. 2007, and references therein). Although there 
are general accepted views on the way star formation occurs 
and young stellar objects (YSO s) evolve to the main sequence 
(IShu et al.lll9"87t lLarsonlll973h . the existing theories have not 
yet converged to an agreed paradigm that can explain the wide 
range of existing observational properties of YSOs. In particu- 
lar, since their discovery, hundreds of brown dwarfs (BD) with 
masses down to the planetary regime have been uncovered in 
star-forming regions and the solar neighbourhood, with a ra- 
tio of the number of BDs to stars of approximately 1/5 (see, 
for example, Luhman et al.l l2007bl and references therein), 
implying that a successful star and planet formation theory 
must account for them. Different theories for the formation of 
BDs are currently debated, according to which they could ei- 



ther form b y gravitationa l fragmentation an d collapse of molec 
ular cores dPadoan et all 120071: iHennebell e & Chabrier 200: ' 
from early ejection from stellar embryos ( Reipurth & Clarki 
2001; Whitworth & Goodwin 2005), o r from f ragmenta tion o f 
massive circumstellar discs (Stamat ellos & Whitworthl l2009l) . 
The extension of the IMF to the brown dwarf and planetary 
mass regime and the search for the end of the mass function 
is therefore crucial to determine the dominant formation pro- 
cess of substellar obje cts and its relat i on with the surround- 
ing environm ent (Mo raux et al.l 120071; lAndersen et alj 120081 ; 
'Luhman"2007t). Brown dwarfs are brighter when they are young 
(Chabrier et al l2000l) and their detection down to a few Jupiter 
masses can be attained with the current techn ology by study 
ing them in y oung sta r-forming r egions (L ucas & Rochd l200 



Zapatero Oso rio et al .112 002; Wei ghts et al.l 2009: 



* Based on observations obtained with WlRCam, a joint project of 
CFHT, Taiwan, Korea, Canada, France, at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institute National des Sciences de I'Univers 
of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii. Based on observations made at the ESO 
La Silla and Paranal Observatory under program 083.C-0092. Based 
in part on data collected at Subaru Telescope, and obtained from the 
SMOKA, which is operated by the Astronomy Data Center, National 
Astronomical Observatory of Japan. Research supported by the Marie 
Curie Research Training Network CONSTELLATION under grant no. 
MRTN-CT- 2006-035890. 



2009MMarsh et al.l2009b . For that reason, one of the prime goals 
of modern observations is to achieve completeness at the lower 
mass end, i.e., the brown dwarf and planetary mass regime, 
for dilfe rent environm ents acros s several you ng star-forming 
regions dBihain et alJ [2 009; Bo uv et alJ l20q9alltt iLodieu et all 



120091; lLuhmanetall2009c iSchoIz et aljr2009l among many oth 
ers). 

The main motivation of our survey of the p Ophiuchi molec- 
ular cloud is to uncover the low-mass population of the cluster 
down to the planetary regime. Despite being one of the youngest 
(~1 Myr) and closest star-forming regions (120 to 145 pc, 
iLombardi et al.l2008l;lMa"maiekl2008h . the high visual extinction 
i n the cloud's core, with Ay up to 50-100 mag ( Wilking & Lad^ 
1983'), make it one of the most challenging environments to 
study low-mass YSOs. The main studies previously conducted in 
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p Oph have been summarised in a recent review (IWilking et alJ 
which includes a census with the ~300 stellar members 
that have been associated with the cloud up to now, from which 
only 15 are estimated to have masses in the substellar regime. 
iMarsh et"al] ([2009i) reported the discovery a young brown dwarf 
with an estimated mass of ~2 -3 Jupiter masses in p Oph, al- 
though we here question its membership to the cloud (see Sect. 
S3. 

We conducted a deep near-IR (7, H, and K^) photometric 
survey centred approximately on the cloud's core and covering 
~ 1 deg^, which we use to identify candidate members in the sub- 
stellar mass regime. Near-IR surveys are particularly suitable to 
study this star-forming region because most of its population is 
visibly obscured. Previous near-IR studies of this cluster have 
been done from the ground down to a sensitivity limit of K < 13- 
14 mag for a larger area of the cloud (Greene & Young| [T992t 
IStrom et all 119951: iBarsonv et all [19971). and of f < 1 5-5 mag 
for a smaller region (200 arcmin . IComeron et al ]|1993h . Deeper 
observations were done from space with a small coverage o f 
72 arcmin^ and a sensitivity of // < 21.5 mag ( Allen et al]|2002h . 
The WIRCam near-IR survey presented takes advantage of a 
new generation of wide-field imagers on 4 meter-class tele- 
scopes, to reach completeness limits of approximately 20.5 in 
J, and 18.9 in // and over the entire degree-size area of the 
sky occupied by the p Ophiuchi central cloud. This work com- 
plements the previous surveys both in the area it covers and in 
sensitivity. Of compar able characteristics is the near-IR survey 
recently conducted by Alves de Oliveira & Casalil (l2008h . which 
uses a different technique, near-IR variability, to select candidate 
members. Our selection method allows BDs with masses down 
to a few Jupiter masses (according to evolutionary models) to be 
detected through ~20 magnitudes of extinction. Extensive use 
of archive data at optical and IR wavelengths is made to further 
characterise the candidate members. In a pilot study, a spectro- 
scopic follow-up of a subsample of these candidates has con- 
firmed six new brown dwarfs. 

In Sects. |2] and [3] the observations and reductions for new 
and archive data are described. Section |4] explains the methods 
used to select candidate members of p Oph and the results, and 
in Sect.|5]we discuss their properties. Section|6]describes the nu- 
merical fitting procedure used to analyse the data from the spec- 
troscopic follow-up and the spectral classification. These results 
are then discussed through Sect. [T] Conclusions are given in 
Sect. El 



2. Observations and data reduction 

We present the deep infrared photometric survey we conducted 
in the p Ophiuchi cluster in the 7, H, and Ks filters, and the in- 
frared spectroscopic follow-up of a small sample of candidate 
members of this star-forming region. 



2.7. The WIRCam/CFHT near-IR survey 

The WIRCam at the CFHT telescope is a wide-field imaging 
camera operating in the near-infrared, consisting of four Hawaii- 
II2-RG 2048 X 20 48 aiTay detectors with a pixel scale of 0.3" 
(iPuget et al.ll2004l) . The four detectors are arranged in a 2 x 2 
pattern, with a total field of view of 20' x 20'. 

The data were obtained in queue-scheduled observing mode 
over several runs as part of a large CFHT key programme aimed 
at the characterisation of the low-mass population of several 
young star-forming regions (RI. J. Bouvier). Seven different 
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Fig. 1. Histogram of the number of objects detected per magni- 
tude bin and the respective magnitude errors. The points where 
the histograms diverge from a linear fit to the logarithmic num- 
ber of objects per magnitude bin give an approximation of the 
completeness limit of the survey for the different filters: 20.5 in 
J, and 18.9 in H and Ks, with errors below ~0. 1 magnitudes. 



WIRCam pointings were needed to cover the central part of the 
p Ophiuchi cluster These were taken at different epochs over a 
three-month period, but all observations were done under pho- 
tometric conditions, with a seeing better than 0'.'8 (measured in 
the images to be typically between 0'.'4 and 0'.'5), and an air- 
mass less than 1.2. All individual tiles were observed in the J, H, 
and Ks filters, using a seven point dithering pattern selected to 
fill the gaps between detectors, and accurately subtract the sky 
background. Table [T] shows the central position (right ascension 
and declination) for each of the seven tiles and the dates of the 
observations. For each field, short and long exposures were ob- 
tained with the J filter (7x4x5 s and 7x8x27 s, respectively), and 
shorter individual exposures with the H and K, filters (7x8x7 s). 

Individual images are primarily processed by the 'I'iwi re- 
ductions pipeline at the CFHT (Albert et al., in prep.), which 
includes detrending (e.g. bias subtraction, flat-fielding, non- 
linearity correction, cross-talk removal), sky subtraction, and 
astrometri c cal ibration. Afterwards, the data are handled by 
Terapix (lMarmoll2007l) . the data reduction centre at the Institut 
d'Astrophysique de Paris (France) responsible for carrying out 
the final quality assessment of the individual images, determin- 
ing precise astrometric and photometric calibrations, and com- 
bining the dither and individual exposures into the final stacked 
images. All images are merged into a single tile of ~1 deg^ 
centred on the cloud's core. The photometric calibration of the 



C. Alves de Oliveira et al.: The low-mass population of the p Ophiuchi molecular cloud 
Table 1. Journal of the WIRCam/CFHT observations. 
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" For the / filter, two sets of images were taken: short (7x4x5 s) and long (7x8x27 s) exposures. 



WIRCam data is done with 2MASS stars in the observed frames 
as part of the nominal pipeline reduction. Typical estimated er- 
rors in the WIRCam zero point determination are of ~0.05 mag. 
Ultimately the photometric accuracy for a determined field is 
also dependent on the number of 2MASS stars available (which 
is probably reduced in regions in the sky with high extinction 
as is the case of some regions of p Oph), if the stars used 
are themselves variable dAlves de Oliveira & Casalill2008l) . and 
lastly, photometric offsets can also occure from small problems 
in flatfielding and/or sky subtraction from night to night, given 
that the data were taken in different epochs. 

We extracted PSF photometry from the mosaicked images 
with PSFEx (PSF Extractor, Berlin et al., in prep.), a software 
tool that computes a PSF model from well-defined stellar pro- 
files in the i mage, which is given as an input to the SExtractor 
programme (iBertin & ArnoutsI 19961) to compute the photometry 
for each detected object. During this first stage of the analysis, 
and to ensure the detection of all the faint sources present in the 
images, the extraction criteria used are not too stringent. An ob- 
ject is extracted if it complies with the required minimum to have 
three contiguous pixels with fluxes 1.5 cr above the estimated 
background. An inspection of the images and detections showed 
that the number of spurious detections is minimum, while all the 
objects seen by eye are detected. Catalogues of the short and long 
exposures for the J filter are merged into one single catalogue. 
The overlap in magnitudes between the short and long exposures 
allows checking of the photometric accuracy. For objects that 
are common to the two catalogues, the two magnitude values for 
each object are compared, and the r.m.s. accuracy is measured to 
be below 0.05 magnitudes. The histogram of the magnitudes (not 
corrected for extinction) is shown in Fig.[T] We derived approxi- 
mate completeness limits of 20.5 in J, and 18.9 in // and Ks, with 
eiTors below ~0.1 magnitudes, located a t the points where the 
histograms diverge f rom the dotted lines (IWainscoat et alJ[T992t 
[Santiago et al.ll 19961) . which represent a linear fit to the logarith- 
mic number of objects per magnitude bin, calculated over the 
intervals of better photometric accuracy. The catalogues from all 
filters are combined into a single database by requiring a posi- 
tional match better than 1" and detections across the J, H, and K, 
lists. The mean separation for the ~27,000 detections common 
to the J band short and long catalogues is found to be ~0'.'05, 
and when combining all the bands is ~0'.' 1 . The final catalogue 
contains ~57,000 objects. Approximately 1000 of the brightest 
stars have a counterpart in the 2MASS catalogues, and the mean 
magnitude differences between the two systems is found to be 
0.05, 0.07, and 0.09 mag for J, H, and Kg respectively, which is 
of the order of expected zero point uncertainties. The dispersion 
of the differences can however be as high as ~0. 1 mag, which re- 



flects the possible sources of error mentioned, and therefore we 
did not correct the photometry for these offsets. Furthermore, the 
WIRCam and 2MASS filters design differ substantially, and at 
no point are these colour effects taken into account. Throughout 
this work, all the WIRCam J, H, and Ks photometry is given in 
the CFHT Vega system. 

2.2. The Soft I NTT spectroscopic follow-up 

A spectroscopic follow-up was conducted for a subsample of 1 3 
candidate members of the p Ophiuchi Cluster, with magnitudes 
ranging from 12.5 to 15 in Ks. The selection of candidate mem- 
bers is described in Sect. |4] We also observed GY 201, which 
was previously associated with the cloud based on its mid-IR 
colours (Wilkin g et al. 200 8) but was not selected with our cri- 
teria. All observations were gathered from 3-6 May 2009, us- 
ing Sofl (Son of ISAAC), a near-IR low resolution spectrograph 
mounted on the 3.6 m New Technology Telescope (NTT, La 
Silla, ESO). The majority of the targets were observed with the 
blue and red grisms, which operate from 0.95 to 1.64, and 1.53 
to 2.52 fj.m, respectively. Some objects were too faint in the J 
band and could only be observed with the red grism. In addi- 
tion, we observed nine sources that are not part of the candidate 
member list, but have magnitudes and colours close to those of 
the selection limits and can therefore serve as a test of our se- 
lection criteria (see also Sect. H)). Field dwarf optical standards 
were also observed (LHS 234, vB 8, vB 10, LHS 2065, Kelu-1), 
though the final spectral classification method we adopted uses 
young optical standards instead (see Sect. |6]l. The observations 
were done with the long slit spectroscopy mode, with a slit width 
of 1" or 2" to better match the seeing conditions, resulting in a 
resolution of ~500 and ~300, respectively, across the spectral 
range. The individual exposure times were chosen according to 
the target's brightness and night conditions, and repeated with an 
ABBA pattern for posterior sky-subtraction. Standard AO stars 
were observed at regular intervals and chosen to have an air- 
mass matching that of the target within 0. 1 . The slit position was 
aligned with the parallactic angle. 

The data reduction was done first with the Sofl pipeline de- 
veloped and maintained by the Pipeline Systems Department at 
the European Southern Observatory (ESO). The 2D spectra were 
flat fielded, aligned, and co-added. The extraction of each spec- 
trum was done with the APALL routine in IRAF. All spectra were 
wavelength calibrated with a neon lamp. The telluric coiTections 
were done by dividing each spectrum by that of the standard AO 
stars observed at similar airmass and interpolated at the target's 
airmass. Relative fl uxes were rec overed with a theoretical spec- 
trum of an AO star (|PickleJl998l taken from the ESO webpage) 
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Fig. 2. Spectra of Kelu-1 (an L dwarf optical standard) taken 
during the observin g run with Soft / NTT (black) and with the 
SpeX / IRTF (grey) dBurgasser et al.ll2007l SpeX Prism Spectral 
Libraries). 



smoothed to the corresponding resolution. To remove the strong 
intrinsic hydrogen absorption lines from the spectra of the AO 
standard stars, a linear interpolation was made across the lines 
that are more predominant at this resolution (the Paschen 6 line 
at 1 .00 /vm, Paschen a at 1 .09 fj.m, Paschen y6 at 1 .28 fim, and the 
Brackett series lines at 1.54, 1.56, 1.57, 1.59, 1.61, 1.64, 1.68, 
1.74, 2.17 fim). The final spectra were not flux-calibrated, but 
simply normalised to their average flux in the 1.67-1.71 fim re- 
gion. The excellent agreement in the overlapping region of the 
spectra taken with the blue and red grisms shows that no further 
calibrations are needed. Figure |2] shows the comparison of the 
SofI / NTT spectrum of Kelu-1 (an L dwarf optical standard) 
taken during this observing run with a sp ectrum of the same 
object taken with the SpeX spectrograph dRavner et alj l2003h 
mounted on the 3 m NASA Infrared Telesc ope Facility, provide d 
in the SpeX Prism Spectral Librarie^ (Burga sser et alJl2007h . 
There is a good agreement in the spectral features and the two 
spectra exhibit similar relative fluxes between the three photo- 
metric bands. 

2.3. The NICS I TNG spectroscopic follow-up 

A shorter observing run was conducted at the Telescopio 
Nacional Galileo (TNG, La Palma, Observatory Roque de Los 
Muchachos) with the Near Infrared Camera and Spectrograph 
(NICS), a low-resolution spectrograph working in the near-IR 
regime. The observations took place during two half nights on 
17-18 May 2009, and three of our candidate members with mag- 
nitudes of Ks~l3 were observed (in addition to the main observ- 
ing programme, which did not directly concern this study), as 
well as the field dwaifs vB 10 and LHS 2924. We used the JK' 
grism, with a wavelength range from 1.15 to 2.23 fim, resulting 
in a resolution power of ~350. The observations were done with 
the 1" slit aligned at the parallactic angle, and AO standard stars 
were observed for telluric corrections. Standard IRAF routines 
were employed to reduce the data, in an analogous way to that 



described in the previous section. The spectra were normalised 
to their average flux in the region between 1.67-1.71 fim. 

3. Archival data 

In order to complement our selection criteria of young stellar ob- 
jects in p Oph and to better characterise the new candidate and 
confirmed members, we made extensive use of multi-wavelength 
data recovered from different archives. The datasets used and 
the criteria employed to extract reliable samples are briefly de- 
scribed in this section. 



3. 1 . Spltzer Space Telescope: C2D surve)^ 

To compleme nt this stud y, Spitzer data from the C2D legacy 
project (E vans et al.ll2003l From Cores to Disks) were included. 
The p Ophiuchi m olecular cloud has been mapped with Spitzer's 
Infrared Camera (i Fazio etal.ll2004 IRAC) in the 3.6, 4.5, 5.8 
and 8.0 /im bands over a region of 8.0 de g^ and with the 
Multiband Imaging Camera (Rie ke et al] l2004. MIP S) in the 24 
and 70 //m bands over a total of 14.0 deg^ tPadgett et al.ll2008l) . 
which encompass the WIRCam field in its totality. The data 
were retrieved from the C2D point-source catalogues of the final 
data delivery (Evans & c2d Team 2005), using the NASA/ IPAC 
Infrared Science ArchiveQ. 

All fluxes were converted to magnitudes using the following 
zero points: 280.9+4.1, 179.7+2.6, 115.0+1.7, 64.1+0.94 (Jy), 
for the 3.6, 4.5, 5.8 and 8.0 //m IRAC bands, respectively, and 
7.17+0.11 (Jy) for the 24 fim MIPS band. Only sources with 
magnitude errors below 0.3 magnitudes as well as detections 
above 2 cr were kept. The Spitzer catalogues were merged with 
the WIRCam detections catalogue, requiring the closest match 
to be within 1". A counterpart was found for ~ 15,000 objects 
that were detected in one or more mid-IR bands. Infrared ex- 
cess around young stellar objects (YSOs) is a direct evidence of 
discs and their detection is commonly used as a youth indicator 
(Haisch et al. 2001). These data are relevant to assess the like- 
lihood of membership for the candidate members as well as to 
characterise their morphological properties (Sect. | 



3.2. Subaru Telescope: i' and z' band archival data 

We searched the Subaru Mita ka Okayama Kiso Archive sys- 
tem (SMOKAlBaba et al.l2002h for Subaru Prime Focus Camera 
(iMivazaki et aL 120021 Suprime-Cam) optical images overlap- 
ping with the WIRCam/CFHT survey. We found one overlap- 
ping field observed in the Sloan /'-band on 20 June 2007, and two 
fields observed in the Sloan z'-band on 16-17 April 2004. Table|2] 
gives a summary of the observations. Each field was observed 
in dithering mode to eff'ectively compute and remove the sky. 
Weather conditions on Mauna Kea on 20 June 2007 and 17 April 
2004 were photometric, as reported b y the CFHT Skyprob e 
atmospheric attenuation measurements (ICuillandre et al.ll2004l) . 
No data are available for 16 April 2004, but the quality and depth 
of the images suggests that the weather was similarly good. 
Seeing as measured in the images ranged from 0'.'6 to ^.'8 in 
the /' and z'-band images. 

The 10 individual CCD of the Suprime-Cam mosaic were 
processed with the standard reduction procedure with the recom- 
mended SDFRED package ( Yagi et al. 2002; Ouchi et al. 2004- 
The programme SDFRED performs overscan and bias sub- 
traction, flatfielding, distortion correction, atmospheric disper- 



Available at |http ://www.browndwarf s . org/spexprism/ 1 



Available at http://irs a.ipac. caltech.edu/| 
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Fig. 3. J-H vs. H-Ks colour-colour diagram (left) and Ks vs H-Ks colour magnitude diagram {right). Plotted over the density map 
for the entire WIRCam catalogue for p Oph are the candidate members selected from near-IR photometry in this study (black filled 
circles), all the candidate members of the cluster from the literature present in the WIRCam catalogues (white circles, asterisks 
signal potential contaminants) including the previously known brown dwarfs from the literature (grey triangles). Candidate members 
selected in this study for the spectroscopic follow-up (squares) and test sources observed spectroscopicall y (crosses) are also s hown. 
In both diagrams, the solid line represents the DUSTY 1 Myr isochrone labelled with solar masses (Mq) (IChabrier et alJl20()Q) . and 
the dashed lines the ~75 and 4 My^^^ limits, with increasing amount of visual extinction. 



sion correction, sky subtraction, masking vignetted regions, and 
alignment and co-addition. A sixth order astrometric solution 
was computed using 2MASS counterparts. The final accuracy is 
expected to be better than 0'.'2. The programme SExtractor was 
used to identify all sources brigther than the 3-cr local standard 
deviation over at least 3 pixels. The absolute zeropoint for each 
CCD was derived from the observation of a SDSS secondary 
standard field (SA 1 10, Schmidt 2002) observed the same night. 
They are given in Table |2] and agree with the iMivazaki et alJ 
(l2002h measurements within 0.12 mag and 0.02 mag in the 
and z'-band, respectively. The chip-to-chip offsets were com- 
puted from the median flux of domeflat images. 



4. Selection of substellar candidate members of 
p Oph 

4. 1 . Colour-colour and colour-magnitude diagrams 

The primary criteria used to select candidate members in the 
substellar regime is to compare the positions of all the WIRCam 
sources in various colour-colour and colour-magnitude diagrams 
with the predictions from the models of the YSOs colours. In 
the first iteration and in an attempt not to exclude any possibly 
interesting candidates, we performed no filtering on the initial 
catalogue of elongated objects (galaxies, nebulosities) or objects 
that could have their photometry affected by instrument artifacts. 



This step was done later with more stringent criteria to ensure 
the quality of the selected candidates. Evolutionary models are 
known to become increasingly uncertain at younger ages and 
lower masses, which is the case for our survey, and a more ac- 
curate way to select candidate members based on photometry 
diagrams is to use o bserved colours of known YSOs instead. 
iLuhman et alJ (1201 Oh empirically determined intrinsic colours 
for young stars and brown dwarfs, which the authors present in 
the 2MASS photometric system. Given the differences between 
2MASS and the WIRCam filter system, which cause large dif- 
ferences in colour, and because to date no colour transformation 
equations between the two systems have been derived, we cannot 
use this approach though. We rely on evolutionary models for 
the selection of candidate substellar members, notin g that when 
comparing the observed colours of YSOs from iLuhman et al.l 
(l20IOh to the Dusty model (2MASS filters, IMyr), we find flie 
differences to be of the order of our photometric colour errors. 
Furthermore, using the models we recover the previously known 
brown dwarfs and the majority of candidate members from the 
literature present in the WIRCam catalogues (see Sect. l4.2l for a 
detailed discussion). 

In the first step, candidate members of p Oph were selected 
if their colours fell redward from the model isochrones in the 
J vs. J-H, J vs. J-Ks, and Ks vs. H-Ks colour-magnitude dia- 
grams. From these, only objects that had colours consistent with 
them being young and substellar in the J-H vs. H-K^ colour- 
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Table 2. Journal of the Suprime-cam / Subaru observations. 



Pointing 


RA 


Dec. 


Date 


Filters 


Exp. Time 


Zeropoint 




(J2000) 


(J2000) 








[ABmag] 


1 


16 27 00 


-24 38 21 


20 June 2007 


Sloan /' 


20x80 s 


28.04+0.03 


2 


16 27 06 


-24 13 30 


17 April 2004 


Sloan z' 


16x240 s 


27.03+0.03 


3 


16 25 06 


-24 38 30 


16 April 2004 


Sloan z' 


21x200 s 


27.04+0.04 



Table 3. WIRCam data for field star contaminants. 



During the spectroscopic follow-up, we observed nine 
.sources that were not part of the substellar candidates list. These 



RA Dec J H K, sources were chosen to pass the selection criteria in the vari- 
(rn^g) (mag) (mag) CMDs, but to have positions in the colour-colour magnitude 



16 25 15.92 
16 25 49.93 
16 26 22.47 
16 26 46.12 



-24 25 10.8 
-24 13 43.7 
-24 37 24.2 
-24 21 53.8 



15.50+0.05 
17.02+0.05 
16.22+0.05 
15.75+0.05 



14.40= 
15.60= 
14.89= 
13.97= 



0.05 
0.05 
0.05 
0.05 



13.46+0.05 diagram near to the reddening line extended from the 75 M /,,„ 
14 73+0 05 " 
13 94+0 05 '^o'l^l colour, which we used for the selection of candidates. The 

12 88+0 05 sources have magnitudes less than ~14 in the Ks band. All these 

~ turned out to be field dwarfs, further supporting the limits used 

in our selection criteria, in particular at this magnitude range. 

These sources are also shown in Fig.^(crosses). 



colour diagram (Fig. [3]l were kept. The 1 Myr Dusty isochrone 
dChabrier et al.ll2000l) computed for WIRCam/CFHT J, H, and 
Kg filters was used down to a temperature of ~1700 K, shifted 
to a distance of 130 pc. The adopted distance to p Oph is a me- 
dian value to that of several distance estimates existing in the 
literature, which indicates that distances to different regions of 
the cloud can v ary from 120 to 145 pc dLombardi et al.ll2008t 
lMamaiekll2008h . According to the models, the substellar limit is 
at /- 11.8, 11.1, and J-H~03. This selection holds 178 sub- 
stellar candidates. From this list, sources were removed which 
had a flux-radius value (flux-radius is a SExtractor parameter 
that measures the radii of the PSF profile at which the flux is 50% 
from its maximum) inconsistent with stellar profiles, i.e. close 
to zero or much larger than the average PSF FWHM measured 
on the images. This ensured that instrumental artifacts (bad pix- 
els, cosmics) and elongated sources (galaxies, nebulosities), re- 
spectively, were correctly discarded. The remaining sources have 
been visually inspected, and a further rejection criterion was im- 
plemented to exclude detections susceptible of having bad pho- 
tometry, like those at the edge of the field, in the overlapping 
regions between detectors (mostly in the detectors edges), or in 
zones of bright reflection nebulae. 

We removed from the candidate list five young brown dwarfs 
(GY 11, 64, 141, 202, CRBR 31, see also Table |7ll, which 
have already been s pectroscopically co nfirmed as members 
(iLuhman et al."1997t (Wilking et al.l[l999t ICushing et al.ll2000t 
[Natta et al. 2002). There are several objects which have previ- 
ously been associated with the cloud but lack a spectroscopic 
confirmation (see also Sect.l4.2|i, in p a rticular from previous IR 
surveys fromlCrreene & Young! (Il992h . IStrom et all (Il995h . and 
iBon temps et aL ( 200 lb: me mbers identified through X-ray emis- 
sion ( Imani shi et al. 200 iL ^Gagne et al. 2004); a nd candidate 
members proposed based on their Spitzer colours dPadgett et alJ 
l200ilWilkingetani2008l:lGutermuth et al.ll2009l see Sect. 15:2]). 
These were kept in our catalogue and are referenced accord- 
ingly when mentioned. We also removed objects from the list 
of candidate members that turned out to be field star contam- 
inants from our spectroscopy follow-up (see Sect. 16.21 the co- 
ordinates and near-IR magnitudes for these sourc es are given 
in Ta ble [3]), as well as three sources observed by iMarsh et"al] 
(l2009i) . which were found not to be substellar (identifiers in that 
study are #1307, #2438, and #2403). The final fist contains 1 10 
substellar candidates selected from near-IR photometry alone, 
which are listed in Table |4] 



4.2. Comparison to previous surveys of fhe p Optiiuciii 
molecular cloud 

To better evaluate the quality of our images, detection meth- 
ods, and the consistency of our candidate selection criteria, we 
compared our results to those of previous surveys of p Oph. 
In a recent compila tion of previous studies in the literature, 
Wil king et al] d2008l) gathered a list of 316 confirmed or can- 
didate members of the cluster, from which 295 have positions 
on sky within the field of our survey. The majority of these are 
part of the initially extracted WIRCam catalogues, with only 13 
objects from the literature missed by the detection algorithm. 
From these, nine sources ( ISO-29, 31, 60, 85, 9 0, 99, 125, 137 
144 (Bontemns et al. 2001); [GY92]-167, 168 ( Greene & YoungI 
1992); HD 147889; CRBR-36 ( Comeron et al. 1993) ) are either 
extremely saturated in the WIRCam images or in the vicinity of 
those and other bright reflection nebulae, and therefore impos- 
sible to detect in the saturated pixels. The remaining four (ISO- 
60, 85, 90, 99) were previously associated with the cloud either 
from their X-ray emission or mid-IR excess. However, there is 
no signal detected in the 7-band of the WIRCam images, hence 
they are not present in the combined JHKs catalogue. Finally, 
sixty-five of the previously known or candidate YSOs are within 
the magnitude range of our near-IR candidate selection, and 35 
of those fall in the substellar region (including the five brown 
dwarfs mentioned in the previous section), further supporting 
our selection criteria. Objects detected by our extraction algo- 
rithm but not within the magnitude range adopted for our can- 
didate selection include seven sources that are too faint in the / 
band and have no reliable magnitude measurement, and 211 ob- 
jects that are brighter than the survey saturation limit for one or 
more near-IR bands. 

While comparing the positions of th e candidate members 
from the literature from the compiled list of lWiUdng et al.l (l2008h 
in the near-IR diagrams, we found inconsistencies between the 
colours of fifteen sources and the expected colours of YSOs. In 
the colour-magnitude diagrams, these sources are in the substel- 
lar regime region (we caution that this is only an approxima- 
tion, given the known errors in the models). However, in the 
colour-colour diagram, they show colours consistent wit h those 
of stars. One of these sources is an edge-on disc (Gros so et al.l 
2003, known as the Flying Saucer). Given the extended pro- 
file of this source, it is likely that its PSF photometry done on 
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the WIRCam images has larger errors, because the parameters 
were optimised for point-source photometry. The remaining 14 
sources include ROX C J162821. 8-245535, which according to 
IWilking et al.l (l2008h has been assigned membership based on 
X-ray emission (unpublished), and 13 sources cla ssified as can- 
didate s from the analysis of IRAC data done by IWilking et alj 
where several diagnostics for the detection of mid-IR 
excess were used (identification numbers used in that work for 
these sources are IRAC 20, 746, 763, 830, 831, 869, 901, 1016, 
1086, 1212, 1343, 1350, 1401). From these, only one source 
(GY 376, or IRAC 746) is present in the l ist of candidate mem- 
bers of p Oph from lGutermuth et alj ('2009) who used the same 
Spitzer dataset, and none ha ve been otherwise pr eviously asso- 
ciated with the cloud. Since IWilldng et alJ dloosh have not pro- 
vided details of their reduction of the data, or mid-IR magnitudes 
for the new candidates we cannot further comment on the valid- 
ity of their selection. However, according to our near-IR dataset, 
these sources seem inconsistent with being members of p Oph 
(see also Sect. 15.21) . 

Additionally, we included in our study recent results not 
found in the com pilation from j^y ilking et al. (2008) from the 
DROXO survey dSciortino et all l2006t Deep Rho Ophiuchi 
XMM-Newton Observation), which consists of a very deep ex- 
posure (total exposure time i s 515 ksec^ taken with the European 
Photon Imaging Camera (Striider et al. 2001; Turner et al. 2001, 
EPIC) on board the XMM-A^ewton satellite ( Jansen et al. 2001), 
and covering a region of ~0.2 deg^ of the p Ophiuchi cluster. The 
data reduction and ma i n results from thi s survey can be found in 
iGiardino et alJ (l2007h . iFlaccomio et al ] 12009), and Pillitteri et 
al. (2010, in prep.). A total of 1 1 1 X-ray emitting sources are re- 
ported in their studies. The sensitivity of the DROXO survey and 
area covered are much lower than that of our survey. We found 
a positional match within 2 cr from the source positional errors 
for two of our candidate members. 

We also c ompared the WIRCam data with the Ks photometry 
presented by iMarsh et"al] (12009) for the seven candidate YSOs 
in p Oph o b served spectroscopic ally in that study (Table |5]). 
iMarsh et"an (l2009 l) derive near-IR photometry from stacking 
deep integration J, H, and Ks images from the 2MASS calibra- 
tion scans. We found a good agreement between their 2MASS 
measurements and the WIRCam photometry for four sources 
(#1449, #1307, #2438, #2403) with magnitude differences be- 
tween 0.02 and 0.23 magnitudes. But for the remaining three 
sources we found larger variations, with differences in mag- 
nitude of 0.4, 1.42, and 1.57 for sources #2974, #4450, #3117. 
Two of these sources (#2974, #31 17), were also detected in 
the WFCAM/UKIRT images from lAlves de Ohveira & CasaUl 
(|2008). and their magnitudes have a difference of 0. 14 and 0. 17 
to those of the WIRCam, which is of the order of their pho- 
tometric errors at this magnitude range and seems to indicate 
that our measurements are correct. Of particular interest is the 
result for source #4450, classified by Marsh et al. (2009) as a 
young T2 dwarf, which may be the youngest and least mas- 
sive T dwarf observed spectroscopically so far. Fro m our images 
we der ived a K, magnitude of 19. 14+0.20, whereas [Marsh et aU 
(l200 9l) reported Ks-ll.Tl. Assuming the parameters derived in 
the spectral analysis of Marsh et al. (2009) are correct and re- 
peating the same calculation as the authors to estimate the dis- 
tance, we arrive at a +lcr range in distance of 137 to 217 pc 
using the WIRCam magnitude value. This seems to suggest 
that this object is behind the p Ophiuchi cloud and could in- 
stead be part of the Upper Sco associa tion ( de Geus et al. 198^ 
located at ~145 pc (de Bruiine et alj jl997l) and an estimated 
age of 5 Myr (IPreibisch & Zinneckei 1999h . As claimed by 




247.2 



246.4 



Fig. 4. Spatial distribution of the substellar candidate members 
of p Ophiuchi in this study and previously known candidate and 
confirmed members of the cloud as compiled by Wilkin g et al.l 
(2008), superposed on the density map of all WIRCam detec- 
tions with contours from the NICER extinction map. Symbols 
are the same as in Fig. [3] 



IMarsh et alJ (l2009l) . in that case the estimated mass for this 
brown dwarf would still be < 3 Jupiter masses, according to 
the models. We could not derive a magnitude for the WIRCam 
H band because the object's position is coincident with an ar- 
tifact caused by the guiding star In the / band, we derived a 
magnitude of 2 1.32 ±0.35. Its J-K, colour is consistent with a T 
dwarf reddened by the extinction measured spectroscopically by 
IMarsh etalJd2009h . 



5. Photometric properties of tlie candidate 
members 

We used multiwavelength data to characterise the candidate 
members presented in Table|4] For each candidate, a flag was in- 
cluded to indicate additional information: previously suggested 
candidate member in the literature (Sect. 14.21 ). candidate selec- 
tion confirmed from optical counterpart (Sect. 15.1b . mid-IR ex- 
cess as determined from Spitzer diagrams (Sect. 15.21 ). variability 
behaviour (Sect. 15.3b . or their membership confirmed spectro- 
scopically from this study (Sect. |6]l. Figure |4] shows the posi- 
tion on sky of the candidate members and the known members 
of p Oph, superposed on the density map of all WIRCam detec- 
tions and the contours from the extinction ma p of the p Ophiuchi 
cloud provided by the COMPLET^ project (Ridge et al."2006'; 
Lombardi et al. 2008), and computed with the NICER algorithm 
dLombardi & AIvesll200lh using 2MASS photometry. 



Available at |http://www.cfa.harvard.edu/COMPLETOl 
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Table 5. WIRCam data for sources in lMarsh et aP (l2009h . 



ID (as in Marsh et al. (2009)) 


RA 


Dec 


J 


H 


a:., 








(mag) 


(mag) 


(mag) 


2438 


16 27 09.37 


-24 32 14.9 


22.21+0.39 


19.53+0.14 


16.97+0.05 


2974 


16 27 16.74 


-24 25 39.0 






17.26+0.06 


3117 


16 27 17.68 


-24 25 53.5 




18.64+0.07 


17.27+0.05 


2403 


16 27 21.63 


-24 32 19.2 


20.86±0.10 


18.17+0.06 


16.50+0.05 


4450 


16 27 25.35 


-24 25 37.5 


21.32±0.35 




19.14+0.20 


1449 


16 27 30.36 


-24 20 52.2 


19.59+0.06 


16.96+0.05 


15.65+0.05 


1307 


16 27 32.89 


-24 28 11.4 


21.41+0.15 


17.44+0.05 


14.97+0.05 




2 4 6 8 10 12 2 4 6 8 10 12 

z-J [mag] i-K [mag] 



Fig. 5. z' vs z'-J {left) and /' vs i'-J {right) c olour magnitude dia grams. In both diagrams the solid line represents the DUSTY 1 Myr 
isochrone labelled with solar masses (M0) (IChabrier et al.ll200Q) . and the dashed lines the ~75 and 4 My^^^ limits, with increasing 
amount of visual extinction. Symbols are the same as in Fig. [3] 



5.1. Near-IR and optical CMDs 

For the candidate members with a counterpart in the /' and z'- 
bands from the Subaru telescope, we could further test our near- 
IR selection criteria. A visual inspection of the match between 
candidate members and the optical counterparts was performed 
to ensure the quality of the positional association. The candidate 
members without an optical counterpart were either to faint in 
the the /' or/and z'-bands or were not in the area covered by the 
optical images. Figure |5] shows the colour-magnitude diagram 
combining infrared and optical data togeth er with the theoret- 
ical isochrone from the DUSTY models OChabrier et alJl2000l) 
for 1 Myr age, shifted to 130 pc. AH but five candidates with ei- 
ther /' and/or z'-band photometry show colours consistent with 
those predicted by evolutionary models within the photometric 
measurement eiTors, and more important with those of the pre- 
viously known brown dwarfs in p Oph. Our selection criteria are 
therefore confirmed for the majority of the sources present in 



these diagrams. Furthermore, GY 201 shows colours that are too 
blue when compared to the isochrones or to the other candidates 
and members. 

5.2. Mid-IR excess from Spitzer data 

Young stellar objects can show infrared emission, which orig- 
inates from dusty envelopes and circumstellar discs suiTound- 
ing the central object. Mid-IR data from the IRAC and MIPS 
Spitzer cameras allow these objects to be studied at wavelengths 
where the excess contribution from discs and envelopes is pre- 
dominant. With several colour-colour and colour-magnitude di- 
agrams (Fig. |6]l we could further characterise our list of candi- 
date members. Previous work on Spitzer mid-IR observations 
of the p Ophiuchi cluster has been published by Padge tt et al.l 
(2008) and Gutermuth et al. (2009), where hundreds of candi- 
date members are uncovered over a much larger area than the 
WIRCam/CFHT survey. 
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Fig. 6. IRAC / Spitzer colour-colour diagrams (panels (a), (b), and (c)), and WIRCam/CFHT and MWSISpitzer colour-magnitude 
diagram. The various diagrams are used to separate young stars of different SED classes, and identify possible extragalactic con- 
taminants. Symbols are the same as in Fig. [3] 



The IRAC colour-colour diagram ([3.6]-[4.5] vs. 
[5.8]-[8.0]) in the panel (a) of Fig. |6] can be used as a 
tool to separate young st ars of different classes jAUen et alj 
120041: iMegeath et alJl2004 and to reject sources consistent with 
galaxies dominated by PAH emission and narrow-line AGN 
(Gutermuth et al. 2009). Centred in the origin are sources which 
have colours consistent with stellar photospheres and have no 
intrinsic IR-excess. These can be foreground and background 
stars, but also Class III stars with no significant circumstellar 
dust. In this region of the colour-colour plane it is impossible 
to differentiate between young stars and contaminants. Another 



preferred regio n for objects in the diagram is located within the 
box defined bv lAUen et alj (l2004i) . which represents the colours 
expected from models of discs around you ng, low-mas s stars . 
Finally, from models of infalling envelopes, lAllen et aP (l2004l) 
predict the colours of Class I sources to have ([3.6]-[4.5]) > 0.8 
and/or ([5.8]-[8.0]) > 1.1. Thirty seven of our candidate 
members have good photometry in the four IRAC bands and are 
displayed in the diagram. The remaining candidates are either 
too faint in the IRAC images or have detections in one of the 
four bands that did not match the quality criteria we applied 
(see Sect. 13. 11 1. From those, 21 were previously associated with 
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p Oph. All the candidates in panel (a) have colours consistent 
with Class 11 or Class 1 young objects. Furthermore, none of the 
candidates falls in the preferred region defined by the colours 
of possible extragalactic contaminants (see Sect. I5.4l i. which is 
also confirmed from the diagram (b). Panel (c) can be used to 
estimate contamination levels by broad-line AGN and is further 
discussed in Sect. 15.41 



Twenty six of the candidate members have a detection at 
24 nm and are displayed in panel (d). Following the lGreene et alj 
( 1199 4^ mid-lR classication scheme (based on the a index), YSOs 
will lay on defined areas of the diagram, namely /T^- [24] > 8.31 
for Class 1, 6.75 < Ks-[24] < 8.31 for flat-specti-um objects, 3.37 
< K,-[24] < 6.75 for Class 11, and K,-[24] < 3.37 for photo- 
spheric colours. All our candidates with a detection at 24 yum 
are consistent with being young according to this classification 
scheme. 



5.3. Near-IR variability of YSOs 

Variability is a characteristic of YSOs, and near-IR variabil- 
ity surveys in particular can probe stellar and circumstel- 
lar environments and provide information about the dynam- 
ics of the on going magnetic and accretion processes. With 
the Wi de Field near-IR camera (W FCAM) at the UKIRT tele- 
scope, lAlves de Ohveira & CasaUl (|2008) conducted a multi- 
epoch, very deep near-IR survey of p Oph to study photometric 
variability. They found 137 variable objects with timescales of 
variation from days to years and amplitude magnitude changes 
from a few tenths to ~3 magnitudes. We found that 17 of our can- 
didates show photometric v ariabihty (14 are in cluded in the list 
of members compiled by Wilki ng et al.l (|2008)), which further 
supports their membership. From their list of candidate members 
found through photometric variability, 18 are outside our sur- 
veyed area, and 58 are saturated in the WlRCam images. From 
the variables that have magnitudes and colours consistent with 
our selection criteria for substellar candidate members of p Oph, 
we recovered all but one. The source AOC J162733.75-242234.9 
has large photometric variations in the near-IR (0.6 and 0.4 mag- 
nitudes in the H and K band, respectively) and is detected in all 
the WlRCam images. But its position overlaps with an artefact 
in the J band image, which may affect its photometry and ex- 
plain the J-H colour found, which is too blue in comparison to 
the theoretical models. From our list of candidates, 93 do not ap- 
pear in the list of near-IR variables. These include 12 candidates 
that a re out of the field surveyed by Alves de Oliveira & Casali 
(l2008h and 24 that are fainter than their completeness limits (~ 19 
and 18 magnitudes, in H and K, respectively). We therefore con- 
clude that the remaining 57 candidate members (from which 13 
have been previously associated with the cloud) did not show 
photometric variability in the near-IR on timescales from days 
to one year during the epochs surveyed in that study. 



5.4. Contamination 

Photometrically selected samples of candidate members of 
young star-forming regions are prone to be contaminated by 
other objects with colours similar to those of YSOs. Only a spec- 
troscopic analysis can reveal the true membership of the can- 
didate members, but this requires large amounts of observing 
time in telescopes and that is not always achievable. The possi- 
ble sources of contamination are extragalactic objects (like AGN 
or PAH galaxies), foreground and background field M dwarfs, 
and background red giants. We tried to estimate the level of 



contamination in our list of candidate members, calling to at- 
tention however that the membership of an individual source is 
ultimately dependent on its position in the field, because for a 
large part of the area of sky surveyed, the cloud's extinction will 
substantially shield any background contamination. We removed 
from the discussion the 30 candidates that were previously asso- 
ciated with the cloud because they gather an ensemble of prop- 
erties from different surveys that further supports their member- 
ship (see Sect. 14.2b . as well as the spectroscopically confirmed 
members in this study, leaving 70 candidates. 

5.4.1. Extragalactic contaminants 

In an attempt to characterize the extragalactic contamination lev- 
els in YSO samples selected on the basis on IR-excess emis- 
sion, i n particular using Spitzer observations, Gutermuth et al.l 
(l2009h explored the same data as lStern et al.l (l2005h to select ac- 
tive galaxies and compare them to YSO selection methods. In 
particular. Stern et al. (2005) found that PAH-emitting galaxies 
have colours that are confined to specific areas in most of the 
IRAC colour-col our diagrams. These regions were adopted by 
iGutermuth et al.l ([2009) to filter out these contaminants, and they 
are depicted in panels (a) and (b) of Fig. |6] None of the candi- 
date members present in this diagram falls into either of these 
regions, indicating that these group of candidates is most likely 
not affected by contamination from star-forming galaxies. 

Another source of contamination are broad-li ne AGN whi ch 
have mid-lR colours similar to those of YSOs ( S tern et al.l2005h . 
With the IRAC diagram panel (c) in Fig. |6] we try to provide 

an estimate for contamination level from AGN, following the 

I 1 ( 

methodology from Guieu et al. (2009) for underredened colours 

(based on Gutermuth et al. 2009). The region plotted in the dia- 
gram (in grey) shows the area in the [4. 5] vs. [4.51-[8.01 colour 
space consistent with AGN-like sources. lGutermuth et al.l (l2009h 
found that while applying this cutoff significantly improved the 
extragalactic filtering of catalogues, some residual contamina- 
tion is still expected. Three of the YSO candidates in these dia- 
gram fall in the contamination area and are signaled out as possi- 
ble contaminants in Table|4] However, only 40 of the candidates 
have detections at 4.5 and 8 pm. If we extrapolate this to the can- 
didate list, we estimate a contamination level of ~5 extragalactic 
sources, a conservative upper-limit because we are not taking 
into account the cloud's extinction. 



5.4.2. Galactic contaminants 

In the Galaxy, giants are an important source of contamina- 
tion. However, taking in account the high galactic latitude of 
the p Oph field (+16.7) and that background contamination is 
reduced due to cloud extinction, fewer giants are expected be- 
cause we are surveying a region of the sky above the plane and 
bulge, which becomes dominated by the faint end of the dwarf 
luminosity function. 

We used the Besan gon model of galactic population synthe- 
sis (i Robin et alJl2003h to estimate the level of contamination by 
foreground late type objects and background red giants or ex- 
tincted galactic sources. We retrieved a synthetic catalogue of 
sources within 0.78^ toward the direction of our survey for dis- 
tances in the range 0-50 kpc. Objects further away than 100 pc 
(hereafter background objects) were placed randomly within the 
field of our survey, and the corresponding extinction as given 
in the COMPLETE map was applied. The luminosities of ob- 
jects closer than 100 pc (hereafter foreground objects) were not 
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Table 6. Spectral type and Ay determined through numerical 
spectral fitting. 



CFHTWIR-Oph 


Spectral Type 
Num. Fit. H2O Index 


Ay (mag) 


4 

34 
47 
57 
62 

96 
106 


M6.50+^;j 
M8.25+f^ 
M7.50+«j 
M7.25+0-^5 

M8.25!*f 
M6.50!|-^^ 


6.41±1 
7.98±1 
7.09±1 
7.12+1 
4.01 + 1 
7.51 + 1 
6.42+1 


^- '"-0.4 
5.6!«i 

6.10:i;' 

1 10+'''' 
4 q+0.6 



changed. Our selection algorithm was then applied to the output 
synthetic catalogue. According to these simulations, only one 
unrelated galactic source could have passed our selection crite- 
rion, indicating that the contamination by galactic sources must 
be low. The contaminant is an extincted (Ay=7.5 mag) 1 Gyrold 
M8V dwarf located at 140 pc. The strong extinction in the region 
covered by our survey indeed most likely blocks the light of the 
majority of background sources up to the limit of sensitivity of 
the survey. 

6. Spectroscopic follow-up of candidate YSOs 

We obtained near-IR spectra for 16 candidate members of p Oph, 
chosen from our WIRCam/CFHT survey, and for GY 201 a can- 
didate member from the literature. The low resolution and mod- 
est signal-to-noise of the spectra (S/N~15, and sometimes lower 
in the J band) restrict the use of narrow spectral features for 
classification. However, even at this resolution, there are still sig- 
nificant differences between the spectra of a low-gravity young 
stellar object and a field dwarf which can be studied. The trian- 
gular shape of the H band, caused by deep H2O absorption on 
either side of the sharp peak located between 1.68 and 1.70 yum 
in young objects, as opposed to a plateau in the spectra of field 
dwarfs, has been used as a signature of youth and member- 
ship in_several studies of young brown dwarfs (see, for exam- 
ple, |AUers_et_aL2007;Lu£as^ There are strong water 
absorption bands on both sides of the peak also in the K band 
spectrum of young brown dwarfs. In the 7-band, H2O absorp- 
tion is also present at both extremes of the band for young brown 
dwarfs, but not for field dwarfs. Another good gravity indicator 
in the near-IR spectrum is the Na I absorption (present at 1.14 
and 2.2 /zm), which is very deep for field dwarfs, but not in young 
objects. Our spectral classification method relies on the compar- 
ison of the candidate spectra with those of young optically clas- 
sified objects members of other star-forming regions of similar 
ages, which are used as standards. A numerical spectral fitting 
procedure was developed, which makes the simultaneous deter- 
mination of spectral type and reddening possible. We also took 
spectra of sources outside the substellar selection limit we im- 
posed in the colour-colour diagram (see Sect. |4]i, to ensure that 
our selection criteria are not too stringent, and indeed all of those 
objects turned out to be stars with no water absorption features. 

6. 1 . Numerical spectral fitting 

The procedure consists in comparing each candidate spectrum 
to a grid of near-IR, low-resolution template spectra of young 
stars and brown dwarfs with spectral types determined in the op- 
tical, reddened in even steps of Ay. The comparison spectra are 



of members of the you ng (<2 Myr ) star-forming regions IC 348 
(^Luhman et al.l l20()3ah and Taurus (iBriceno et al.ll2002L iLuhmanI 
[2004) dereddened by their extinction published values (typically 
Ay < 1), with spectral types ranging from M4 to M9.5. By com- 
bining the spectral types, half a nd quarter sub-clas ses were con- 
structed. The reddening law of iFitzpatrickl (1 19991) was used to 
progressively redden the template spectra by steps of 0. 1 Ay. 
The fit was performed across the complete usable wavelength 
range of the spectrum (1-1.34, 1.48-1.8, and 2-2.45 jum), i.e., 
excluding only the regions dominated by telluric absorption and 
the extremes of the spectral range where the quality of the data 
is poorer It was assumed that all template spectra have approx- 
imately the same error For the candidate spectra, the r.m.s. of 
the difference between the original spectrum and the smoothed 
spectrum (using a 15 pix boxcar) was taken as an estimation of 
the errors. Figure [T] shows the contours for the variation of 
with Ay and spectral type for one of the candidates. The solid 
line contour represents the 1 sigma confidence interval, and the 
dashed lines indicate the Ay and spectral type limiting values at 
the point of the grid closest to the 1 sigma contour, which are 
taken as the standard deviation for each parameter from the min- 
imum. The dotted contours are the 1.6, 2, 2.6, and 3 sigma levels, 
successively from the minimum. The right-hand panel shows the 
resulting best fit, CFHTWIR-Oph 34 is best-fitted by an M8.25 
brown dwarf (an average spectra between an M8 and M9 young 
brown dwarfs) and an Ay of 9.7 magnitudes. This procedure was 
applied to all the candidate spectra. 

6.2. Spectral classification: comments on individual sources 

6.2.1. Contaminant field stars 

Water vapour absorption could not be detected in the spectra of 
four candidates, and only a limit to the spectral type could be 
set, i.e., they have a spectral type earlier than ~M-type. Given 
their faint IR magnitudes, they are therefore inconsistent with 
being young and members of the cluster and are excluded as 
background contaminants (listed in Table[3]l. 

6.2.2. CFHTWIR-Opli 4, 34, 47, 57, 62, 96, 106 

The sources CFHTWIR-Oph 4, 34, 47, 57, 62, 96, and 106 have 
spectral types and extinction values derived from the numerical 
procedure, with spectral types ranging from M5.5 to M8.25, and 
Ay from a 1 to 10 magnitudes. The results are summarized in 
Table|6] while the dereddened spectra are displayed in Fig.|8] 

6.2.3. GY 201 

The source GY 201 was observed with the blue and red grisms 
of Sofl/NTT, and could neither be fitted with one of the tem- 
plates in the grid, nor with comparison spectra from field dwarfs 
when the full spectral range was considered. If the fit was per- 
formed with only the part of the spectrum acquired with the red 
grism (from 1.5 to 2.5 yum), the fitting procedure converged for 
spectral type M5 and an Ay of 1.9 magnitudes. However, when 
the full spectrum was taken into account, no physical solution 
was found, with the resulting best fit indicating a negative value 
of Ay. This disagreement can be explained if the object is an 
unresolved binary, where one of the components is an earlier 
type dwarf contributing to the blue part of the spectrum, and 
the other a late M-type dwarf dominating the red part of the 
spectrum and showing water vapour absorption features. Yet an- 
other plausible and more likely explanation is that the red part 
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Fig.7.x^ map obtained in the spectral type vs. Ay plane for one of the candidates (CFHTWIR-Oph 34), together with the best-fit 
solution. The spectrum shows a very good match to that of an intermediate spectrum between an M8 and M9 young brown dwarf 
of Taurus (IMyr). 



of its spectrum is contaminated by the emission from a nearby 
(<10") Class I young member of p Oph previously studied by 
many authors (ISO 103, see for example, Gutermuth et al.'2009^, 
[Padgett et al. 2008; Imanishi et al. 2001; Bontemps et al. 2001). 
The source ISO 103 is very bright in the K band, and is satu- 
rated in our WIRCam images. The source GY 201 was also pre- 
viously associated with the cloud based o n two near-lR studies 
(iGreene & Young|ll992l:lAilen et al.ll200l . but was not detected 
or mentioned in any other study of the cloud, even if its location 
has been covered by the vast majority of the surveys. Its position 
on the various optical and near-IR CMDs presented in this pa- 
per indicates that this object could be a field dwarf, because its 
position is always bluewards of the models and other candidate 
members. It also lacks an IRAC or MIPS detection in the Spitzer 
archive catalogues, which could be explained by the difficulty of 
separating its PSF from that of the neighbouring Class I source, 
which is also very bright in the m id -IR. It has also not been de- 
tected by Gutermuth et al. (2009) or Padgett et al. (2008) in their 
two detailed Spitzer studies of p Oph. The nature of this object, 
and in particular its association to the cl oud, remains therefor e 
uncertain. Furthermore, in the review bv 'Wilking et al.1 (l2008h . 
the authors claim this object to be a Class I source based on (un- 
published) IRAC colours, a result which seems unlikely taking 
into account the results from the spectrum analysis in our work 
and its near-IR colours, and could be a mismatch between the 
IRAC detections and the existing literature catalogues, or IRAC 
photometry confusion caused by the neighbouring star 

6.2.4. CFHTWIR-Oph 2, 55, 94, 97, 105 

The candidate members CFHTWIR-Oph 2 and 94 show a very 
red spectrum, without photospheric features required for a spec- 
tral classification, like clear water vapour absorption bands. The 
nature of these sources is further discussed in Sect. [7] Figure [8] 
shows the original spectra not dereddened, because they lack a 
classification. These objects were observed with a grism cover- 
ing the entire near-IR spectrum (NICS/TNG), while CFHTWIR- 



Oph 55, 97, and 105 were observed with Sofl/NTT using only 
the red grism, because they are too faint in the J band. These 
spectra are also very red, though water vapour absorption can 
be seen in the spectra of CFHTWIR-Oph 97 and 105. These 
candidate members are classified as M5.5 and M6, respectively, 
though their classification is less reliable given that only a lim- 
ited part of their spectra is available for the fit. The two other 
objects did not show clear absorption bands and have not been 
fitted. Their original undereddened spectra are also displayed in 
Fig- El and their properties are discussed in detail in Sect.|7] 

6.3. The H2O Spectral Index 

For the seven candidate members with a spectral classifica- 
tion we compared the results from the numerical spectral fit- 
ting process with the H2O spectral index defined bv .AUers et al.l 
(2007), which can be calculated as (/^,i=i.55o-i.56o/-F^=i.492-i.502)- 
Allers et al. (2007) derive a spectral type vs. index relationship, 
which is independent of gravity and is valid for spectral types 
from M5 to LO, with an uncertainty of + 1 subtype. The index 
is computed for all our spectra after being dereddened by the 
Ay values in Table|6l and the spectral types all agree with those 
determined by the fitting procedure within the uncertainties and 
are shown for comparison, further confirming the validity of our 
classification method. 

7. Properties of the spectroscopic sample 

7.1. Membership 

We used the compiled information for each candidate member 
observed spectroscopically to confirm their pre-main-sequence 
nature. For the seven sources that have determined spectral types 
and extinction values (Table |6l), we found an agreement with the 
spectra of young stellar objects with low gravity, and therefore 
identify in their spectra signatures of youth like the //-band trian- 
gular shape. Additionally, all objects that have an optical coun- 



C. Alves de Oliveira et al.: The low-mass population of the p Ophiuchi molecular cloud 13 




nl I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 



Fig. 8. Sofl/NTT and NICS/TNG low-resolution spectra of the observed candidate (left panel) and confirmed (right panel) members 
in p Oph. The spectra on the right panel are corrected for extinction with the values found through the numerical fitting procedure, 
which are shown together with the spectral type derived for the best-fit solution. 
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terpart (CFHTWlR-Oph 34, 47, 62, 96, and 106) show colours 
similar to those of the previously known brown dwarfs in p Oph. 
In the mid-lR diagrams, CFHTWIR-Oph 2, 34, 55, 62, 94, 96, 
97, and 105 show evidence of discs (Sect. lS!2] ). 

Some of these objects have been previously associated 
with the cloud. The source CFHTWIR-Oph 62 was previ- 
ously associated with the cloud from a comparison of near- 
IR photometric observations to model s and the detec t ion of 
infrared excess ("Riek e & Riek3 1 19901: IComeronetaP Il993t 
[Greene & Young 1992), and it was first observed spectroscop- 
ically by Wilking et al. ( 199^ but lacked a high enough signal- 
to-noise to be studied. Cus hing et al.l (l2000 l) observed the same 
object in the near-lR (NIRC/Keck I), and claim its member- 
ship based on the detection of strong H2 emission in the K- 
band, though the authors mention it is not clear if the emission 
is associated with the object. H2 emission is also detected in 
our Sofl/NTT spectra, but the resolution of the spectrum is too 
low for an accurate velocity measurement ( s ee als o Sect. 17.21 ). 
The spectral type derived bv ICushing et alj (|2000) is M4±1.3, 
which agrees within the errors, with the spectral type we found, 
M5.5+'^^. The sources CFHTWIR-Oph 34 and 96 have also 
been previously associated with the cloud based on the detection 
of near-lR excess (Greene & Young 1992), but lack a spectro- 
scopic confirmation. These three sources have mid-IR colours 
consistent with those expected for YSOs dWilking et al J 120081 
iGutermuth etani2009h . Finally, CF HTWIR-Oph 34 and 62 have 
been c lassified as variable sources bv Alves de Ofiveira & CasaUl 
(l2008h . which further supports their classification as members. 
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Fig. 9. SED of a candidate edge-on disc in p Oph. The SED of 
CFHTWIR-Oph 62 is compared to the SED of a known Class 111 
young stellar object in p Oph with a similar spectral type, red- 
dened by an extinction amount to match that of the candidate. 
Both sources are scaled at the //-band flux. This object is un- 
derluminous at optical wavelengths, does not show an excess 
emission from the near-lR up to 8 //m, but has a large excess at 
24 jdm. 



Despite the fact that the candidate members with very red 
spectra could not be classified, there is evidence they are mem- 
bers of the cluster. The sources CFHTWIR-Oph 55, 94, 97, and 
105 have been pre viously associated w ith the cloud from IR and 
/ or X-ray surveys dWilking et al.ll2008h . The source CFHTWIR- 
Oph 2 is a new candidate member and shows colours consistent 
with those of a Class 11 object (Fig.|6]i. 



7.2. H2 Outflow 



The source CFHTWIR-Oph 94 (other names are, for example, 
GY 312 or ISO 165) is a known member o f p Oph and has been 
extensively studied: llmanishi et al.' (2001') detected both quies- 
cent and flare X-ray emission, Natta et al. (2006) fo und it to 
be an actively accreting YSO, and Alves de Oliveira & Casalil 
(l2008l) detected photometric variab ility consistent with ch anges 
in the surrounding disc or envelope. lBontemps et alJ (1200 ih clas- 
sified it as a Class 11 object, i.e., with an IR excess and a spec- 
tral energy distribution (SED) which can be explained by mod- 
els of YSOs surro unded by circumstellar discs. More recently, 
using Spitzer data, IGutermuth et al.l (l2009l) have classified it as 
a Class 1, given its strong IR excess. Our mid-IR colour-colour 
magnitude diagrams agree with the later classification. We found 
further evidence of the protostellar nature and therefore youth of 
this object. We detected a H2 1 - S(0) emission (2.12 /^m) in the 
spectrum of CFHTWIR-Oph 94 (Fig. [8]), a signature of a molec- 
ular outflow. Given the extreme red spectrum of this object, we 
cannot estimate its spectral type. Further observations are needed 
to determine the association of the outflow with the source (see , 
for example, iBourke et al.ll2005t [Fernandez & Comer6nll2005h . 
and also its mass. 



7.3. Candidate edge-on disc 

We investigated the mid-IR colours of a candidate edge-on disc, 
CFHTWIR-Oph 62, which shows very red colours at 24 //m 
but is not present in the IRAC diagrams (the object is detected 
at 3.6, 4.5, and 5.8 pm, but has only an upper limit detection 
at 8 yum). We compared the colours of CFHTWIR-Oph 62 to 
those of WSB 50, a young stellar object member of p Oph with 
a spec tral type close to th at of CFH TWIR-O ph 62, between 
M4.5 dWilking et alJ l2005h and M4 dLuhman & Rieke 1999h . 
and classified as a Class III source, which should therefore show 
a photospheric SED. The magnitudes are reddened to those of 
CFHTWIR-Oph 62 (with an Ay of 8., WSB 50 has an Ay^l.9 
magni tude s) with the reddening laws from iRieke & LebofskyI 
dl985h and lFlahertvetal] d2007l) . The near-IR magnitudes were 
taken from the 2MASS catalogue because WSB 50 is saturated 
in the CFHT/WIRCam images, and it is preferred to use a com- 
mon photometric system. Figure |9] shows the two SEDs, which 
nicely shows the ch aracteristic SED of CFHTWIR-Oph 62. 
ICushing et alJ d2000l) found this object to be underluminous at 
optical wavelengths, which combined with the lack of IR ex- 
cess up to 8 pva could be explained by the geometry of a nearly 
edge-on disc: at short wavelengths the disc is optically thick 
and acts as a natural coronograph (explaining the underlumi- 
nosity when compar ed to members of similar Tf/y, and why 
ICushing et alJ (l2000h concluded this object to be older), while at 
longer wavelengths the thermal emission of the disc domi nates, 
causi n g the sharp rise in the SED (see, for example Saut er et al.l 
120091: (Puchene etani2009l) . That the flux at 24 pm is at the 
level of the photospheric part of the SED further supports this 
scenario. Further observations and modelling are needed to un- 
derstand and better characterise this complex young object. In 
particular, measurements at longer wavelengths than 24 yum can 
provide an important constraint in the nature of this source. 
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Fig. 10. H-R diagram for p Oph displaying the new brown dwarf candidates found i n this study, togeth er with previously known 
stellar and substellar members of the cloud. The models are the NextGen isochrones (iBaraffe et al.lll998h shifted to the distance of 
the cluster (130 pc), for 1, 5, 10, and 30 Myr, and labelled with mass in units of Mq. 



7.4. Temperatures and luminosities 

To place the new candidate members that were spectroscopi- 
cally confirmed in the Hertzsprung-Russell (HR) diagram, which 
is commonly used to estimate ages and masses, we needed 
to derive their effective temperature and bolometric luminos- 
ity. To convert spectral types to temperatures, the temperature 
scale from J^uhman et al. (2003b) was adopted, which is de- 
rived for young members of the star-forming region IC 348 
(~2 Myr), and has provided consistent results w hen applied to 
other young star-forming regions (for example, iLuhman et al.l 
12009 ). The adopted eiTors for the temperature are the one sigma 
limits in spectral type from the numerical fitting procedure. 
For each candidate, the bolometric luminosity was calculated 
from the dereddened J magnitude (using A y derived from the 
numer ical fit and the reddening law from Rieke & Lebofskvl 
(Il985l) ). applying t he bolometric correcti o ns for the respective 
spectral type (fro m iKenvon & Hartmaniil (Il995h for <M6 and 
iDahn et alJ (l2002h for >M6), and using a distance to the cloud 
of 130 pc. The errors were propagated to include the photo- 
metric error in J, the one sigma errors in Ay from the spec- 
tral fitting, and by assuming an error of +10 pc in the distance 
to the cloud. For comparison purposes, we compiled the pre- 
viously known members of p Oph with assigned spectral types 
later then > M3 from the literature. From the -300 objects asso- 
ciated with the cloud, members that have a counterpart and are 
not saturated in the WIRCam/CFHT 7-band images were kept. 
Some of these objects are not part of the final WIRCam cata- 



logue (see Sect. 121), because they are saturated either in the H or 
the Ks bands. Only objects that have spectral types determined 
from spectroscopic surveys and extinction values published are 
included. The final co mpilation contains 36 young low-mass 
stars and brown dwarfs (Luhman et al."1997'; 'Luhman & Riekd 
1999; Wilkins et al. 1999; Gushing et al. 2000; Natta et al .l2002l; 
Wilking et al. 2005). Bolometric luminosities and temperatures 
were derived in the same way as for the GFHTWIR-Oph candi- 
dates. The previously confirmed members and candidate mem- 
bers from this study were placed in the HR diagram (Fig.fTOli and 
compared to theoretical evolutionary models (NextGen, because 
our candidates have T^^ > 2500 g. lBaraff^e et al.lll998l) . The 1, 
5, 10, and 30 Myr isochrones are shown, labelled with mass in 
units of Mq. We adopted the 0.08 Mq mass track as the stellar / 
substellar boundary that corresponds to spectral types ~M6.25 to 
M6.5 for a young member of p Oph with an age of 1 to 2 Myr. In 
our sample, we find that GFHTWIR-Oph 62 is a very low-mass 
star, and the other sources are six new brown dwarfs of p Oph 
(GFHTWIR-Oph 4, 34, 47, 57, 96, 106). 

A large spread in ages is seen in the HR diagram, which 
goes from <1 to ~10 Myr, and with some objects laying al- 
ready closer to the 30 Myr isochrone . The estimated age for 
p Oph is of 0.3 My r in the core jGreene & Meveii 119951; 
Luhman & Rieke 1999), and 1-5 Myr in the suiTounding regions 
(Bouvier & Appenzeller 1992; Martin et al. 1998; Wilking et al| 
2005). The star formation history of this cluster is thought to be 
connected to that of the Sco-Gen OB association, with two dif- 
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ferent episodes of star formation taking place, one caused by a 
supernova 1 to 1 .5 Myrs ago, and the other happening in parallel 
to the formation of Upper Scorpius ~5 Myrs ago, caused by an 
exp anding shell from the Upper Centaurus-Lupus OB subgroup 
(see IWilking et alj l2008 l and references therein for a review). 
These ideas are still debated, but would mean that a range of ages 
is therefore expected in the HR diagram. Taking into account 
the typical uncertainties involved in the temperature and Ay de- 
termination, that would explain the position of the members of 
the cluster from above the 1 Myr isochrone up to the 10 Myr 
models. Most of the members in the HR diagram are consistent 
with this age estimate, which is also the case for our candidates 
CFHTWIR-Oph 34, 57, and 96. The other CFHTWIR-Oph can- 
didates have luminosities that suggest an age older than 10 Myr 
up to 30 Myr. We already mentioned, however, that CFHTWIR- 
Oph 62 is underluminous (Gushing et al. 2000), most probably 
due to observations being done through scattered light from a 
surrounding disc (Sect. 17.31 ). The old ages implied in the di- 
agram for CFHTWIR-Oph 4, 47, and 106, do not seem plau- 
sible if we assume them to be members of the cluster. The 
source CFHTWIR-Oph 106 shows mid-lR colours consistent 
with those of young objects surrounded by discs, and it is possi- 
ble that it is observed through scattered light, which would ex- 
plain its lower luminosity. Both CFHTWIR-Oph 4 and 47 do 
not show a signature of IR excess and could be more evolved 
young objects. As mentioned in Sect. 16.11 their spectra are well 
fitted by those of young objects, showing distinctive features of 
youth, in particular the triangular shape of the //-band. FigurefTTI 
shows the dereddened spectrum of CFHTWIR-Oph 4 together 
with the best-fit spectrum from the young grid of templates, 
and a comparison spectrum of a field dwarf of the same spec- 
tral type, putting into evidence the pronounced and broad wa- 
ter a bsorption bands associated with low surface gravity objects 
dKirkpatrick et al. 2006). The same check was done for all our 
classified spectra. Furthermore, three other brown dwarfs taken 
from the literature (CRBR 3 1,GY 11, GY 141) fall into the same 
part of the HR diagram. We do not find any relation between the 
different positions of brown dwarfs in the HR diagram (younger 
or older than the 10 Myr isochrone) and their positions on sky 
in relation to the cluster's core, for example. Nor do we find 
a relation between their position in the HR diagram and their 
SED class assigned from mid-IR colours. In particular, all but 
one of the previously known brown dwarfs in the cloud with an 
assigned SED class are Class II objects. 

The old ages implied by the isochrones could instead be re- 
lated to the several sources of error associated with this diagram, 
like the less reliable photometry for sources located in high neb- 
ulosity regions (abundant in p Oph) , near-IR variability in YSOs 
(lAlves de Oliveira & Casalil (l2008l) detected near-IR variations 
as large as 0.5 magnitudes for objects plotted in the HR dia- 
gram), unresolved binaries, the large difference in estimated dis- 
tances to the cloud, or the possibility that some objects may be 
seen through scattered light. Although the individual contribu- 
tion of these uncertainties can be quantified, it is not possible to 
have a clear picture of the net effect when one or more of the 
mentioned problems are involved. F urthe rmore, recent results in 
modelling of young brown dwarfs (iBaraffe et alJI2009l) suggest 
that episodic strong accretion might explain the observed spread 
in HR diagrams at ages of a few Myr years, a scenario supported 
by recent observations of protostars, some of which were car- 
ried out in p Oph (Enoch etal. 2009). According to these re- 
sults, even after accretion has halted, young low-mass objets can 
keep a memory from these strong accretion events, altering the 
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Fig. 11. Dereddened spectrum of CFHTWIR-Oph 4 together 
with the best fit obtained. The spectrum shows a very good match 
to that of an intermediate spectrum between an M6 and M7 tem- 
plate young brown dwarfs (solid line) and clear differences to 
the spectrum of the fi eld dwarf 2MASS J 1327239 1+0946446 
(iBur^asser et al.ll2004t) with a comparable spectral type (dotted 
Une). 

expected path in their contraction along the Hayashi track, and 
therefore their position in the HR diagram. 

Another possible explanation could also be that some of 
these brown dwarfs lie behind p Oph and are instead members 
of Upper Sco, which would mean the luminosity is underesti- 
mated in the HR diagram presented. If we compute the bolo- 
metric luminosities using a distance of 165 pc instead (approxi- 
mate boundary of Upper Sco), all sources are within the 10 Myr 
isochrone or younger in the HR diagram. Though it is unlikely 
this is the case for all the sources, it is possible that some of the 
brown dwarfs associated with p Oph could rather be Upper Sco 
members. 

Given these uncertainties, we have therefore not assigned a 
definite age or mass to the newly confirmed brown dwarfs and 
the very low-mass star discovered in this work. We can claim 
though, based on their position in the HR diagram relative to the 
other members of the cloud, that they have ages and luminosities 
which agree with those of the known substellar members. From 
their location in the HR diagram, these new candidates indeed 
appear to be amongst the lowest mass objects of the cluster. 

7.4.1 . p Ophiuchi: census update of the substellar population 

We compiled from previous studies a list of spectroscopically 
confirmed members of p Oph with spectral types later than ~M6 
that are therefore likel y to be brown dwarfs (according to the 
evolutionary models of iBaraffe et al.l (Il998h ) and present them 
in Table [T] together with the six new brown dwarfs found in this 
work. All surveys were conducted in the main cloud, L1688. 
All but one (not in the WIRCam/CFHT survey coverage) of 
the brown dwarfs are plotted in the HR diagram in Fig. [TOl 
We did not i nclude in this list the brown dwarfs discovered by 
lAllers et aP (l2007h in a region to the north west of the central 
cloud, because it has been suggested that two o f the three are 
associated with an older population of Sco-Cen dLuhman et al.l 
l2007ah and therefore their membership to the p Ophiuchi cloud 
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Name 


J 

(mag) 


Ay 
(mag) 


Sp. Type 


(K) 


Kei. 


CFHTWIR-Oph 4 


14.88 


2.5 


M6.50 


2935!J| 
2795+^^: 


this work 


CFHTWIR-Oph 34 


15.97 


9.70 


M8.25 


this work 


CFHTWIR-Oph 47 


15.94 


5.6 


M7.50 


this work 


CFHTWIR-Oph 57 


15.06 


6.10 


M7.25 


2880:'™ 


this work 


CFHTWIR-Oph 96 


14.60 


1.10 


M8.25 


27 r 0+310 


this work 


CFHTWIR-Oph 106 


15.36 


4.9 


M6.50 
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5 


oph-193 


13.61'' 
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" PSF photometry magnitude from WIRCam/CFHT survey, unless noted otherwise. 
* From the SpT- l^jj relation of lLuhman eTakl j2003bh . 

Spectral typ es determined from the follo wing studies: 1 . iLuhman et al.) ( Il997h : 2. IWilking et alj il999l) : 3. iLuhman & Riekd (Il999l) : 
4. ICushing et "SI ( 1200 0); 5. Natta et al. ( 2003) ; 6. IWilking et^ ( 120051) . 

This object is not part of the WIRCam surveyed region. Photometry is from 2MASS public catalogue. 



complex has not been confirmed (see also, IClose e t al?'2007'). 
We di d not include in this list the T2 dwarf found by|Marsh et al. 

either, because the main argument used to claim mem- 
bership to p Oph relies on the distance determination using the 
apparent A'-band magnitude of the object, which we claim might 
be wrong. This list provides an updated census of the substellar 
members of p Oph known to date. 

8. Conclusion 

We identify 110 substellar candidate members of p Ophiuchi 
from a deep, near-IR photometric survey, from which 80 were 
not previously associated with the cloud. By extensive use of 
archive multi-wavelength data, we find evidence of mid-IR ex- 
cess for 27% of the candidates and a variability behaviour con- 
sistent with that of YSOs for 15%, further supporting the mem- 
bership of these candidates. 

We started a spectroscopic follow-up of the substellar candi- 
date members, and present the first results for 16 sources. We 
identify six new members of p Ophiuchi with spectral types 
ranging from ~M6.5 to ~M8.25, and classify them as new con- 
firmed brown dwarfs according to the evolutionary models of 
Baraffe et al. (1998). We confirm the spectral type derived by 
Gushing et al. (2000) for a previously known very low-mass star 
close to the substellar limit, and based on the SED constructed 
from optical to mid-IR photometry, we report the discovery of 
a candidate edge-on disc around this star. We cannot derive ac- 
curate spectral types for five sources which have extremely red 
spectra. Two of these show water absorption features and are 
classified with spectral types M5 and M6. However, since they 
lack a y-band spectra and given the poor fit they remain as can- 
didate members. The remaining three sources could be T Tauri 



star members of the cluster, because they show strong mid-IR ex- 
cess and one of them is emitting in X-rays. We found signatures 
of outflow activity in two of the sources studied spectroscopi- 
cally where H2 1 - S(0) emission (2.12 //m) was detected. Four 
sources out of the 16 were found to be contaminant field dwarfs. 
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Table 4. Candidate Members of the p Ophiuchi Molecular Cloud 
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13.34±0.05 
17.25±0.05 
15.12±0.05 
15.68±0.05 
13.11±0.05 
15.32±0.05 
16.69±0.05 
17.31±0.05 
15.25±0.05 
16.13±0.05 
16.4 ±0.05 
13.11±0.05 
13.55±0.05 
14.61±0.05 
16.04±0.05 
15.86±0.05 
15.8 ±0.05 
13.83±0.05 
14.54±0.05 
17. 11 ±0.05 
17.45±0.05 
17.39±0.05 
17.38±0.05 
15.94±0.05 
16.36±0.05 
13.46±0.05 
13.06±0.05 
12.66±0.05 
17. 11 ±0.05 
15.68±0.05 
13.48±0.05 
16.52±0.05 
15.69±0.05 
13.18±0.05 
16.28±0.05 
17.07±0.05 
13.22±0.05 
14.76±0.05 
16.08±0.05 
12.66±0.05 
12.92±0.05 
17.27±0.05 
13.45±0.05 
13.97±0.05 
15.8 ±0.05 
13.3 ±0.05 
14.04±0.05 
12.87±0.05 
14.39±0.05 
14.92±0.05 



M6.5o:0:g 2.5-«:| 



15.01±0.08 14.81±0.09 

11.7 ±0.06 10.96±0.06 
14.06±0.06 13.9 ±0.06 13.83±0.17 
15.2 ±0.08 14.83±0.09 
15.73±0.11 15.28±0.12 

15.76±0.14 15.29±0.1 



9.64 ±0.1 6.88±0.16 



14.35±0.1 



13.93±0.07 13.44±0.( 



.06 


13.41±0.17 




.06 


11.99±0.06 


10.89±0.06 


.05 


10.25±0.06 


9.69 ±0.1 


.06 


11.23±0.08 


10.66±0.07 


.06 


11.34±0.06 


10.75±0.06 7.75±0.14 


.09 






.07 






.06 


11.13±0.24 




.07 


14.24±0.29 




.08 
.05 


10.66±0.06 


10.0 ±0.06 6.62±0.11 


.06 


13.58±0.13 


13.28±0.27 


.09 






.05 


10.88±0.06 


10.66±0.24 


.06 


12.15±0.08 





1,8 
3,8 



M8.25!0;g 9.70!°;J 3,5,8 



1,3,4,7,8 



3,7,8 
4,8 



10.23±0, 
11.33±0, 
15.81±0, 
10.81±0, 
13.35±0, 
14.07±0, 
10.56±0, 
10.87±0, 
11.13±0, 
14.36±0, 
13.61±0, 



06 9.45 ±0.05 8.82 ±0.05 5.07±0.13 

05 10.95±0.06 10.6 ±0.06 10.58±0.06 
08 15.28±0.09 

06 9.84 ±0.06 9.25 ±0.05 8.53 ±0.05 5.6 ±0.1 

07 13.17±0.07 13.05±0.11 
13.28±0.1 



).73 ±0.05 9.29 ±0.05 
10.0 ±0.06 9.34 ±0.05 



12.96±0.2 
13.33±0.28 
8.99 ±0.05 
8.82 ±0.05 
10.0 ±0.05 9.34 ±0.06 



06 13.65±0.06 
06 
06 

06 10.51 ±0.05 

07 14.31 ±0.08 

06 13.14±0.06 12.59±0.08 11.94±0.07 9.22±0.25 



M7.50+°J 



5 6+''-3 
•'•"-0.2 



4.89d 
6.08d 
5.85d 



:0.1 
:0.1 
:0.1 



6,8 



1,3,8 



1,3,4,! 
1,3,4,! 
1,3,4,! 

4,8 



exl,ex2 
exl,ex2,R 



exl,ex2 
ex2 



var,exl,ex2 



exl 
var,exl 
var,exl 
exl,ex2 



exl,ex2 
exl,AGN? 



exl 



var,ex2 



var,exl,ex2 

exl,ex2 
var,exl,ex2 

exl,ex2 



n 



O 



i 



o 



var,ex2 AOC J162624.29-241§9.7 



s 



O 

"a 



2MASS J16263682-24§002 
BKLT J162638-243W1 



[SKS95] 162338.6-24g51 



BBRCG 1 



s 

p. 

o 



BKLTJ162642-2424a9 
BKLTJ162643-242^2 

GDS J162648.3-242^g4 



var,exl,ex2 2MASS J16265197-24M394 



2MASS J16265346-2432362 
BKLTJ162654-242440 
[GY92] 154 



Table 4. continued. 



CFHTWIR-Oph 


/' 


z' 


J 


H 


Ks 


[3.6] 


[4.8] 


[5.4] 


[8.0] 


[24] 


SpT 


Ay Literature 


Flag 


Other name 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(mag) 







J162658.47. 
J162658.66. 
J162703.59- 
J162704.02- 
J162705.94- 
J162707.70- 
J162709.03- 
J162709.80- 
J162710.03- 
J162710.17- 
J162713.04- 
J162713.24- 
J162714.34- 
J162715.69- 
J162715.88- 
J162718.53- 
J162719.39- 
J162719.45- 
J162719.67- 
J162721.ll- 
J162722.41- 
J162724.10- 
J162724.39- 
J162725.64- 
J162726.23- 
J162726.61- 
J162727.65- 
J162730.55- 
J162730.69- 
J162731.74- 
J162732.13- 
J162732.55- 
J162732.73- 
J162734.13- 
J162734.47- 
J162736.04- 
J162736.59- 
J162737.03- 
J162737.21- 
J162737.22- 
J162738.95- 
J162740.12- 
J162740.84- 
J162741.61- 
J162744.20- 
J162745.77- 

100 J162746.54- 

101 J162747.25- 

102 J162750.77- 

103 J162810.46- 

104 J162811.61- 

105 J162821.70- 

106 J162829.93- 



54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 



242004, 
■242455, 
■242005, 
240246, 
241840, 
■243403, 
243025, 
■243442, 
242913, 
■243545, 
243200, 
■242347, 
■243132, 
■243845, 
■242514, 
■240722, 
242600, 
242049, 
■244148, 
■243753, 
■243837, 
242510, 
■244147, 
■243728, 
241923, 
■244045, 
■243827, 
241456, 
■244417, 
■243148, 
■242943, 
241604, 
244500, 
243308, 
241439, 
■244325, 
245136, 
■244334, 
■243434, 
■242526, 
244020, 
242636, 
242900, 
■244644, 
■235852, 
■244453 
240559, 
■244645, 
■244245, 
■242420, 
■243729, 
■244247, 
245406, 



5 
6 

6 22.83±0.02 20.41±0.03 
9 

3 20.62±0.03 
5 
3 


4 20.71±0.01 19.48±0.03 
9 
4 
4 

22.64±0.03 
7 
2 

4 


9 18.68±0.03 
8 
9 
9 
8 

.6 22.22±0.03 

1 23.56±0.03 20.46±0.03 
3 
2 
5 
9 
7 
8 
5 
5 
9 
7 
4 

1 19.85±0.03 
9 
2 
7 

7 22.27±0.0120.13±0.03 
6 

8 18.96±0.01 



18.7 ±0.03 
20.32±0.03 

22.08±0.03 
21.51±0.03 



23.49±0.05 
21.86±0.03 



20.55±0.08 
20.22±0.07 
17.13±0.05 
15.06±0.05 
17.05±0.05 
22.33±0.33 
22.14±0.38 
22.44±0.31 
16.55±0.05 
22.01±0.21 
22.13±0.39 
20.88±0.12 
18.42±0.05 
18.76±0.05 
20.23±0.08 
14.76±0.05 
21.71±0.34 
21.73±0.21 
16.02±0.05 
22.83±0.44 
22.2 ±0.27 
22.32±0.7 
18.78±0.05 
18.23±0.05 
16.35±0.05 
18.49±0.05 
21.94±0.21 
21.75±0.21 
21.57±0.18 
21.92±0.24 
18.35±0.05 
22.42±0.33 
19.03±0.05 
22.22±0.26 
21.93±0.71 
22.84±0.49 
16.83±0.05 
21.91±0.22 
21.92±0.14 
22.31±0.36 
16.48±0.05 
21.9 ±0.19 
14.6 ±0.05 
17.87±0.05 
17.08±0.05 
17.74±0.05 
17.94±0.05 
19.26±0.05 
22.87±0.41 
17.74±0.05 
19.51±0.05 
16.66±0.05 
15.36±0.05 



18.12±0.06 
17.15±0.05 
15.01±0.05 
13.92±0.05 
15.77±0.05 
19.0 ±0.08 
19.3 ±0.12 
18.84±0.07 
15.14±0.05 
18.92±0.07 
19.28±0.11 
17.6 ±0.05 
16.53±0.05 
15.22±0.05 
15.9 ±0.05 
13.55±0.05 
18.71±0.08 
18.35±0.06 
15.15±0.05 
19.47±0.09 
19.19±0.08 
18.73±0.07 
15.01±0.05 
16.54±0.05 
14.36±0.05 
14.89±0.05 
18.94±0.07 
18.73±0.07 
18.82±0.07 
18.47±0.06 
15.04±0.05 
18.81±0.07 
15.41±0.05 
18.96±0.07 
19.05±0.08 
19.05±0.08 
15.65±0.05 
19.16±0.08 
18.21±0.05 
18.81±0.07 
14.12±0.05 
17.39±0.05 
13.76±0.05 
15.87±0.05 
15.78±0.05 
14.8 ±0.05 
16.34±0.05 
16.48±0.05 
19.42±0.08 
16.15±0.05 
16.95±0.05 
14.15±0.05 
14.51±0.05 



16.68±0.05 
14.84±0.05 
13.72±0.05 
13.18±0.05 
14.85±0.05 
16.58±0.05 
16.15±0.05 
16.51±0.05 
14.23±0.05 
17.1 ±0.05 
17.05±0.05 
15.67±0.05 
15.3 ±0.05 
12.95±0.05 
13.26±0.05 
12.79±0.05 
16.79±0.05 
16.17±0.05 
14.57±0.05 
17.36±0.05 
16.92±0.05 
15.85±0.05 
12.72±0.05 
15.36±0.05 
13.04±0.05 
12.54±0.05 
16.95±0.05 
16.92±0.05 
17.06±0.05 
16.44±0.05 
13.05±0.05 
16.57±0.05 
13.25±0.05 
17.03±0.05 
17.24±0.05 
16.71±0.05 
14.85±0.05 
17.0 ±0.05 
15.97±0.05 
16.18±0.05 
12.56±0.05 
14.13±0.05 
13.19±0.05 
14.22±0.05 
14.98±0.05 
12.61±0.05 
15.26±0.05 
14.78±0.05 
17.18±0.05 
15.07±0.05 
15.41±0.05 
12.48±0.05 
13.85±0.05 



15.64±0.07 15.29±0.11 

12.2 ±0.06 11.35±0.05 10.77±0.06 10.1 ±0.06 8.03±0.15 
12.11±0.06 11.45±0.05 10.83±0.06 10.08±0.06 5.91±0.1 
12.53±0.07 12.37±0.07 12.14±0.08 11.96±0.09 



13.75±0.06 13.27±0.06 



12.03±0.11 



M7.50+°| 



5 6+0-3 

-'•"-0.2 



13.6 ±0.06 12.95±0.06 12.44±0.07 12.02±0.13 
14.84±0.06 14.41±0.07 14.25±0.19 
13.16±0.06 12.73±0.06 12.26±0.11 



6.85±0.11 M5.50+0|5 9.90+i| 



15.3 ±0.06 
14.18±0.06 
14.07±0.06 

11.03±0.06 
12.3 ±0.07 
15.32±0.08 
14.68±0.06 
13.65±0.08 
15.72±0.09 
15.22±0.07 
13.77±0.06 
10.93±0.06 
14.16±0.06 
11.64±0.06 

15.3 ±0.07 
15.56±0.08 
15.45±0.07 
15.05±0.06 
11.33±0.06 
15.03±0.07 
11.63±0.06 
15.53±0.07 
15.84±0.08 
15.1 ±0.07 
13.86±0.06 
15.62±0.08 
14.17±0.06 
14.43±0.06 
10.62±0.06 
11.62±0.06 
12.51±0.06 
11.32±0.05 

9.7 ±0.06 
14.41±0.07 
13.35±0.06 
15.6 ±0.07 
13.98±0.06 
13.81±0.06 
10.84±0.05 
13.05±0.06 



14.93±0.07 
13.84±0.06 
13.49±0.06 
9.45 ±0.08 
10.22±0.05 
12.04±0.06 
14.87±0.08 
14.26±0.06 
13.6 ±0.12 
15.44±0.11 
14.91±0.08 
13.22±0.06 
10.47±0.06 
13.9 ±0.06 
11.14±0.05 

14.87±0.08 
15.35±0.09 
14.94±0.07 
14.66±0.07 
10.73±0.05 
14.68±0.08 
11.18±0.06 
15.23±0.07 
15.58±0.13 
14.72±0.07 
13.53±0.05 
15.26±0.11 
13.75±0.06 
13.97±0.06 
9.57 ±0.05 
10.93±0.05 
12.24±0.06 
10.28±0.05 
13.32±0.07 
8.71 ±0.06 
14.17±0.08 
13.01±0.06 
15.32±0.08 
13.68±0.06 
13.22±0.06 
10.24±0.05 
12.57±0.05 



14.33±0.2 
13.44±0.14 
13.06±0.1 

9.6 ±0.05 
11.86±0.07 



11.9 ±0.09 7.2 ±0.11 

6.05 ±0.09 

8.78 ±0.05 5.4 ±0.1 

11.84±0.1 



14.29±0.28 13.57±0.25 
13.2 ±0.28 

14.53±0.29 
12.8 ±0.1 

10.0 ±0.06 10.11±0.06 
13.52±0.11 13.42±0.19 
10.62±0.06 9.93 ±0.05 7.79±0.12 

14.79±0.28 



15.01±0.28 
14.57±0.3 
10.25±0.06 
14.21±0.25 
10.91±0.05 
14.53±0.24 

14.3 ±0.18 
13.15±0.08 

13.39±0.11 
13.73±0.13 

8.7 ±0.05 
10.38±0.05 
11.79±0.06 
9.38 ±0.05 

12.98±0.1 
7.97 ±0.05 
13.96±0.15 
12.72±0.07 
14.95±0.27 

13.29±0.1 
12.62±0.07 

9.7 ±0.05 
12.11±0.07 



9.69 ±0.06 6.22±0.12 
10.79±0.06 8.9 ±0.17 

12.32±0.07 8.91±0.28 

13.66±0.23 

7.7 ±0.05 2.98±0.1 
10.03±0.06 7.67±0.12 
10.96±0.06 8.33±0.17 M8.25 
8.38 ±0.05 5.04±0.1 
12.28±0.11 
7.27 ±0.06 3.88±0.1 

12.65±0.11 



12.97±0.12 
11.92±0.07 
9.23 ±0.05 6.18±0.1 
11.47±0.07 



+0.25 
1.0 



1.10 



+0.7 
0.1 



1,8 
4,8 

2,3,4,8 



var,exl,ex2,R AOC J162658.65-242455.5 
var,exl,ex2 ISO-Oph94 

n 



exl 
var,ex2 



[WGM99] 2408.6-2: 



var,exl,ex2,AGN? AOC J162714.34-243&1.9 
WL20S H 
[GY92] 241 g" 



1,3,4,7,8 var,exl,ex2 



exl,AGN? 



1,2,3,5,8 

1,4,8 
3,5,8 



5,8 
5.8 



1,3,4,8 
4,5,8 
3,4,5,8 
1,8 

1,3,4,8 



M6.50+!-^5 4.9+°-? 



var,exl,ex2 

exl,ex2 
ex2 

exl,ex2 



var,ex2,R 
exl,ex2 
exl,ex2 
var,exl,ex2 

ex2 



exl 
exl 
exl,ex2 
exl 



[AMD2002] JI627I9-2I26OI 



o 

CRBR 2422.6-350g- 
a 

ISO-Ophl38 ^ 
2MASS J16272661-24#451 

O 

"a 



BKLT J162732-242»3 
3 

BKLT J162732-244^ 

s 
s 

p. 

o 

s 

2MASS J16273894-24'|)206 
GDS J162648.3-242§4 
BKLTJ162740-242Sgl 
BKLTJ162740-242^1 

BKLTJ162745-244'§4 



BKLTJ162821-244246 



Table 4. continued. 



CFHTWIR-Oph i' 


z' J H K, [3.6] [4.8] [5.4] [8.0] [24] SpT 


Ay Literature 


Flag 


Other name 


(mag) 


(mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) 


(mag) 






107 J162848.71-242631.8 


14.31±0.05 13.6 ±0.05 13.13±0.05 12.55±0.()6 12.26±0.06 11.88±0.06 1 1.29±0.06 9.73±0.28 




exl,ex2 


BKLT J162848-242631 


108 J162857.87-244055.1 


18.74±0.05 15.93±0.05 14.01±0.()5 11.15±0.()6 10.06±0.05 9.13 ±0.05 8.06 ±0.05 4.12±0.1 




exl,ex2 


SSTc2d J162857.9-244055 


109 J 162902.95-244040. 8 


21.64±0.17 18.98±0.08 17.08±0.05 16.03±0.09 15.7 ±0.11 








110 J162903.97-244105.1 


22.76±0.45 19.38±0.1 17.39±0.05 16.01±0.09 15.51±0.09 






n 
> 



Q 

'<' 

p 

i 
3 

I 

S 

o 

O 

I 
& 

3 
o 

I 

n 
s 

o 

s 

re, 

I 



2 



J2000.0 lAU designation. a_ 
Aperture photometry magnitudes from PrimeCam / Subaru. 
PSF photometry magnitudes from WIRCam / CFHT. 

IRAC and IVIIFS / Spitzer data as retrieved from the NASA/ ffAC Infrared Science Archive. 
^ Spectral Type and Ay as determined from this study. See Sect. |6.1| for details. 

^ p Ophiuchi members , ac cording to the follovy ing studies: l.lBontemns et al.l ( 120011) : 2. IComeron et aLl( fT993h : 3. iGreene & Youria ( fl992l) : 4. iGutermuth etai] l l2009h : 5. iStrom et al.] ( fT995h : 
6. IWilking&Ladal ( ll983h : 7. IWilking etai] i ll 9891) : B. IWilking et al.l ( l2008h . 
^ exl, ex2 = mid-IR excess; Var = Variability; AGN? = Possible Active Galactic Nuclei Contaminant; R = red spectrum. See Sects. [5] and |6.1| for details. 
^ WL 20S, component of the triple system WL 20. 



